Thermodynamic Properties of Heisenberg Ferrimagnetic Spin Chains: 
Ferromagnetic- Antiferromagnetic Crossover 
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We study thermodynamic properties of the one-dimensional Heisenberg ferrimagnet with antifer- 
romagnetically exchange-coupled two kinds of spins 1 and 1/2. The specific heat and the magnetic 
susceptibility are calculated employing a modified spin-wave theory as well as a quantum Monte 
Carlo method. The specific heat is in proportion to T 1 / 2 at low enough temperatures but shows a 
Schottky-like peak at mid temperatures. The susceptibility diverges as T~ 2 . We reveal that at low 
temperatures the model is regarded as a ferromagnet, while at mid temperatures it behaves like a 
gapped antiferromagnet. 

PACS numbers: 75.10.Jm, 65.50.+m, 75.40. Mg, 75.30.Ds 



A great progress has been made in studying the qual- 
itative difference jlj between the integer-spin and the 
half-odd-integer-spin Heisenberg antiferromagnets. Re- 
cently there has appeared brand-new attempts to ex- 
plore the quantum behavior of mixed-spin chains with 
two kinds of spins. Low-energy properties of various 
mixed-spin chains with singlet ground states were ana- 
lyzed Q via the nonlinear-cr-model technique with par- 
ticular emphasis on the competition between the mas- 
sive and the massless phases. Mixed-spin chains with 
magnetic ground states serve us with another topic and 
have in fact attracted much current interests |5j-|| . Since 
we expect a gapless excitation from the ferrimagnetic 
ground state, we there take little interest in the most 
naive problem whether the spectrum is gapped or gap- 
less. Performing a numerical investigation by the use 
of conformal invariance, Alcaraz and Malvezzi H indeed 
predicted the appearance of quadratic dispersion rela- 
tions for the mixed-spin Heisenberg ferrimagnets, which 
is consistent with a spin- wave calculation |i|,|5| . Actually 
the quadratic dispersion has explicitly been visualized 
employing a quantum Monte Carlo (QMC) technique || 
and an exact-diagonalization method [g. Thus we may 
expect quantum ferrimagnets to behave like ferromagnets 
at low temperatures. However, several authors |^,|| have 
recently reported that quantum ferrimagnets have a non- 
trivial excitation branch gapped from the ground states 
as well as a gapless one, which stimulates us to investi- 
gate their thermodynamic properties. We show in this 
article that the two distinct low-lying excitations result 
in a novel temperature dependence of the thermal quanti- 
ties displaying both ferromagnetic and antiferromagnetic 
aspects. The present study sounds more fascinating con- 
sidering that all the mixed-spin-chain compounds syn- 
thesized so far exhibit ferrimagnetic ground states || . 

We consider alternatively aligned two kinds of spins 
S and s on a ring with antiferromagnetic exchange cou- 



pling between nearest neighbors, which are described by 
the Hamiltonian 



N 



(1) 



where M = S z + s z = J2j=i( S j + s j) is the total mag- 
netization, N the number of the unit cells, \ib the Bohr 
magneton, and we have set the g factors of the spins 
S and s equal to g. For the sake of argument, we as- 
sume throughout the manuscript that S > s. Apply- 
ing the Lieb-Mattis theorem [ L0| to the Hamiltonian (jl]) 
with no external field, we find(5 — s)A r -fold degenerate 
ground states. Therefore the model exhibits ferrimag- 
netism instead of antifcrromagnetism. The gapless and 
the gapped excitations, respectively, lie in the subspaces 
of M < (S - s)N and M > (S — s)N and thus may be 
regarded as ferromagnetic and antiferromagnetic. In the 
case of (S, s) = (1, 1/2), the gap to the antiferromagnetic 
branch was numerically estimated to be 1.759 J (111 . 

Now we present QMC calculation of the thermal quan- 
tities at zero field. The recent field-theoretical argument 
H and density-matrix renormalization-group (DMRG) 
study Q both suggest that the low-temperature proper- 
ties of the model are qualitatively the same regardless of 
the values of S and s as long as they differ from each 
other. That is why we restrict our numerical investiga- 
tion to the case of (S, s) = (1, 1/2). We employ the QMC 
method based on the Suzuki- Trotter decomposition |Tl| l 
of checkerboard type and its numerical procedure 
has been detailed elsewhere [Q. We mainly calculate 
the N — 32 chain, which is long enough to discuss the 
bulk properties. Since the correlation length of the sys- 
tem is smaller than the length of the unit cell ^||| , the 
thermal quantities show no significant size dependence. 

We show in Fig. [I temperature dependences of the 
specific heat C (a) and the magnetic susceptibility \ (b) ■ 
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Although the overall temperature dependences of the re- 
cent DMRG findings are similar to ones of our QMC 
results, the two calculations are not in quantitative agree- 
ment with each other. We have confirmed that quantum 
transfer-matrix calculation for short chains precisely 
reproduce the present QMC findings except for very low 
temperatures, where overshort chains may pretend to be 
gapped. Furthermore high-temperature series-expansion 
calculation helps us to verify our numerical treatment. 
Within the up-to-i -3 approximation, the specific heat 
and the magnetic susceptibility are expanded as 



C 

xJ 

,2, ,2 



= f J + 0(i~ 4 ) . 



11 



11 



12 



* + 36 r + °( f ) > 



(2) 
(3) 



where t = k B T/J with the Boltzmann constant k B . The 
asymptotic curves (||) and (||) are also shown in Fig. [l], 
which convincingly fit the numerical results. 

The usual spin- wave treatment diagonalizes the Hamil- 
tonian (@) with no field as |,| H = E s + J2k(^k a l a k + 
w kPkftk), where u>^ = ujk T J(S — s) with uj k = J[(S — 
s) 2 + AS s sin 2 (ak)] 3 , E g = E Q + E\ is the ground state en- 
ergy with E = -2JSsNJ and E x = - J(S + s)], 
and a is the lattice spacing. a\ and p\ are the creation 
operators of the ferromagnetic and the antiferromagnetic 
spin waves with momentum k. The spin-S* ferromagnetic 
Heisenberg chain exhibits the spin-wave excitations with 
a quadratic dispersion u) k — 2JS [1 — cos(afc)] . Thus, 
only in the S — 2s cases, the ferromagnetic branch of the 
spin- (5, s) ferrimagnets show exactly the same dispersion 
as the spin-s ferromagnets exhibit at small momenta in 
the unit of the unit-cell length being unity. Hence we 
expect the present model to behave like the spin- 1/2 fer- 
romagnet at low temperatures. 

The precise low-temperature behavior of the spin-1/2 
fcrromagnet has been revealed by Takahashi and Yamada 
Jl5|. Numerically solving the thermodynamic Bethe- 
ansatz integral equations, they succeeded in expanding 
the thermal quantities by powers of t 1 ! 2 as 



C 
Nfo 



0.7815t3 - 2.00i + 3.5*3 + o{t 2 ) . 



(4) 



x / „ = o.04i67r 2 + o.i45r* + o.m -1 + 0(f*), (5) 

Ng 2 ^ 

which are also plotted in Fig. |l|. Although the 
QMC calculation can not reach low enough tempera- 
tures, yet our findings allow us to conclude that the 
present model is identified with the spin-1/2 ferromagnet 
at low enough temperatures. We note that the lowest- 
temperature QMC estimates, which successfully imply 
the T 1 / 2 asymptotic behavior of the specific heat and the 
T~ 2 divergence of the magnetic susceptibility, were ob- 
tained through the improved algorithm jl3| by spending 
forty million MC steps on each data point. 



At mid temperatures in the specific heat, the antiferro- 
magnetic aspect most clearly appears. The specific heat 
exhibits a sharp peak, rather than a broad one charac- 
teristic of ferromagnets, at k B T j J ~ 0.74 and therefore 
reminds us of the Schottky anomaly peculiar to the anti- 
ferromagnetic specific heat Jl^,|lJ,|l^, ItJ . It is interesting 
to fit the QMC result to the Schottky-type specific heat 



C 
Nk Y 



= A 



A 



2k B T 



sech 



A 



2k B T 



(6) 



with A being set equal to the excitation gap to the anti- 
ferromagnetic branch of the present model, 1.759 J |y§, 
and a fitting parameter A. We find a fine fit with A = 1.7 
as shown in Fig. [j] and thus recognize that the mid- 
temperature behavior of the specific heat is well attributed 
to the gapped antiferromagnetic excitations. 

Now the ferromagnetic and the antiferromagnetic as- 
pects of the model are both revealed. We inquire fur- 
ther into this picture developing the modified spin-wave 
(MSW) theory plf . Introducing the additional con- 
straint of the total magnetization being zero into the the- 
ory, Takahashi jl9| not only overcame the difficulty in the 
conventional spin- wave theory but also succeeded in cor- 
rectly evaluating the low-temperature behavior of various 
thermal quantities. His idea was further applied to the 
two-dimensional isotropic antiferromagnets [ p0pl| ] and 
thus opened the way for a quantitative argument of low- 
dimensional magnets in terms of the spin-wave picture. 
At finite temperatures, we replace atcKfe and /3t/3fe in the 
spin-wave Hamiltonian by rv^ = ^2 n - n + n ± Pfc(n _ ,n + ), 
where Pk(n~ ,n + ) is the probability of n~ ferromagnetic 
and n + antiferromagnetic spin waves appearing in the 
fc-momentum state and satisfies J2 n - n+ Pk(n~~ ,n + ) = 1 
for all fc's. Then the free energy at zero field is expressed 
as 

k 

+ fc BT^ Pk(n-,n+)lnP k (n-,n+). (7) 

k n~ ,n+ 

First, we consider minimization of the free energy (j?]) 
with respect to Pk(n~ , n + )'s under the condition of zero 
magnetization, 

(S z + s*) = (S- s)N - - n+) = . (8) 

k 

The free energy and the magnetic susceptibility at 
thermal equilibrium are obtained within a set of self- 
consistent equations: 

F = E g + M (5 - s)N -kBT^Y, Ml + "*) - (9) 



k a=± 
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3k B T 



(10) 



k <y=± 
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with = [e^i- ± f J -)/ k B T _ l]- 1 ^ where fx is a Lagrange 
multiplier determined by the condition (||). The suscep- 
tibility has been obtained by calculating the thermal av- 



erage of M 2 pi 
in powers of tP 



. Equations (^|) and ( |To| ) are expanded 
at low temperatures and result in 
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f x + 0(f), 



(11) 



(12) 



where C(z) is Riemann's zeta function. In the case of 
(£, s) = (1,1/2), the expressions ( pi] ) and (12) coincide 
with Eqs. (||) and (||) up to the order t and the order t _1 , 
respectively, and therefore again show us £/ie identity be- 
tween the present model and the spin-1/2 ferromagnet at 
low temperatures. 

The above-demonstrated MSW approach gives a sat- 
isfactory description at low temperatures, whereas we 
immediately find that it never applies away from the 
low-temperature region. For ferromagnets, the zero- 
magnetization constraint is not only convincing in that 
the thermal average of the magnetization should be zero 
at zero field, but also plays a role of keeping the number 
of bosons constant. This is not the case for ferrimagnets 
as well as for antiferromagnets pO|,^l| . Here the spin- 
wave treatment with the condition (g) still results in the 
divergence of the number of bosons at high temperatures. 
Can we control the number of bosons keeping the above- 
obtained low-temperature behavior unchanged and give 
a convincing description in a wider temperature region? 
We may answer yes replacing the condition (||) by 



(: S z 



:) = (S + s )N-J(S + s)J2Y,^= - 



(13) 



where the normal ordering is taken with respect to a 
and [3. The spin-wave theory shows us that the clas- 
sical staggered magnetization (S + s)N is modified into 
(S+s)N — 2t with a quantum spin reduction r pj. Equa- 
tion ( |l3| ) claims that the thermal fluctuation of the stag- 
gered magnetization be constrained to take the classical 
value. This is analogous to Eq.(||), which claims that the 
thermal fluctuation of the magnetization be the classical 
magnetization (S — s)N. We stress that the constraint 
(13) leads in fact to exactly the same expressions as Eqs. 
(11) and ( p"2| ) at low temperatures. Now we again obtain 
a set of self-consistent equations: 

F = E g + »(S + s)N -k B TY,J2 ln (! + "*) • ( 14 ) 



k <J=± 



3k B T 



(15) 



k <r=± 



with n 



( u ± Uk -^(S+ s ))Mfe B T _ where // is 



Lagrange multiplier due to the condition ( jjji| ) 



We numerically solve Eqs. (|13| ) and (|14|) in the ther- 
modynamic limit, and visualize them for (5, s) — (1, 1/2) 
in Fig. ||, where the QMC estimates are shown again. 
We find that the MSW calculation not only correctly de- 
scribes the actual behaviors at low enough temperatures 
but also well reproduces the overall temperature depen- 
dences. The T 1 / 2 standing up, the Schottky-like peak, 
the T~ 2 decay of the specific heat, and the T~ 2 diver- 
gence, the T -1 decay of the susceptibility, they are all 
successfully described by our MSW approach. However, 
the spin-wave excitations underestimate the peak tem- 
perature of the specific heat. This is because the spin- 
wave theory results in the gap (jJ^ =0 = J, which is smaller 
than the true value 1.759 J. In order to separately observe 
the contributions of the ferromagnetic and the antiferro- 
magnetic modes, we regard Faf = ~kBTJ2k m (l + "fc ) 

and xaf = (sMb^^bT) -1 Y^k™t (1 + ™k) as tne anti " 
ferromagnetic contribution, while we define the ferromag- 
netic background as Fp = F — Faf and xf = X ~ Xaf- 
Each contribution in the specific heat is numerically cal- 
culated from Fp and Faf , respectively. We show in Fig. 
H the thus-obtained separate contributions. We clearly 
observe that Cf — * C and \f — ► X as T — > 0. On 
the other hand, due to the excitation gap, Caf and xaf 
exponentially vanish as T — > 0. It is due to the anti- 
ferromagnetic mode that the Schottky-like peak appears. 
Finally, while the above consideration is enlightening, we 
admit that the present definition for Cf and Caf may 
not be relevant at high temperatures, where the /x-term 
in Eq. (|l4j) should not simply be incorporated into the 
ferromagnetic part. Such a ambiguity inevitably occurs 
because here the thermal quantities are obtained through 
the nonlinear equations. 

We have investigated thermodynamic properties of the 
fcrrimagnetic mixed-spin chains and revealed that the 
ferromagnetic and the antiferromagnetic aspects simulta- 
neously lie in the model. It may also be emphasized that 
the spin-wave theory has so successfully been applied to 
the model of one dimension as to describe its whole ther- 
mal behavior. Not only direct observation of the thermal 
quantities but also neutron-scattering measurements of 
the antiferromagnetic excitations are fascinating experi- 
ments worth trying. 
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FIG. 1. Quantum Monte Carlo Calculation (O) of the 
specific heat (a) and the magnetic susceptibility (b) as a 
function of temperature for the Heisenberg ferrimagnetic 
spin chain with alternating spins 1 and 1/2 of iV = 32. 
The solid lines: High-temperature series-expansion results 
for the present model within the up-to-( J/k^T) 3 approxima- 
tion. The dashed lines: Low-temperature Behavior [15] of the 
spin- 1/2 ferromagnetic Heisenberg chain obtained by solving 
the thermodynamic Bethe-ansatz integral equations, where 
the thermal quantities are expanded by powers of {k-gT / J) 1 ^ 2 ■ 
The specific heat and the magnetic susceptibility are plotted 
within the up-to-(fcB?7 J) 3,/2 and the up-to-( J/kgT) approx- 
imations, respectively. The dotted line: The Schottky-type 
specific heat with its gap being set equal to 1.759J, that is, 
the energy difference between the ground state and the anti- 
ferromagnetic excitation branch of the present model. 



FIG. 2. Modified spin-wave calculation (solid lines) of the 
specific heat (a) and the magnetic susceptibility (b) as a 
function of temperature for the Heisenberg ferrimagnetic spin 
chain with alternating spins 1 and 1/2 in the thermodynamic 
limit, where both ferromagnetic and antiferromagnetic spin 
waves contribute to the thermal quantities. We effectively 
extract each contribution of the ferromagnetic and the antifer- 
romagnetic spin waves and show them by dashed and dotted 
fines, respectively. Quantum Monte Carlo calculation (O) is 
again plotted for the sake of comparison. 
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Fig. 1 (a) 
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